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We present a detailed analysis of the effective force between two smooth spherical colloids floating 
at a fluid interface due to deformations of the interface. The results hold in general and are applicable 
independently of the source of the deformation provided the capillary deformations are small so that 
a superposition approximation for the deformations is valid. We conclude that an effective long- 
ranged attraction is possible if the net force on the system does not vanish. Otherwise, the interaction 
is short-ranged and cannot be computed reliably based on the superposition approximation. As an 
application, we consider the case of like-charged, smooth nanoparticles and electrostatically induced 
capillary deformation. The resulting long-ranged capillary attraction can be easily tuned by a 
relatively small external electrostatic field, but it cannot explain recent experimental observations 
of attraction if these experimental systems were indeed isolated. 

I. INTRODUCTION 

In view of various applications such as the study of two-dimensional melting [1] , investigations of mesoscale structure 
formation [2] or engineering of colloidal crystals on spherical surfaces [3] , the self-assembly of sub-^m colloidal particles 
at water-air or water-oil interfaces has gained much interest in recent years. These particles are trapped at the 
interface if the colloid is only partially wetted by both the water and the oil. This configuration is stable against 
thermal fluctuations and appears to be the global equilibrium state, as it is observed experimentally that the colloids 
immersed in the bulk phases are attracted to the interface [1] (see also Sec. II). The mutual interaction between 
the trapped colloids at distances close to contact, i.e., within the range of molecular forces, is dominated by strong 
van-der-Waals attraction. In order to avoid coagulation due to this attraction, the colloids can be stabilized sterically 
with polymers or with charges such that the colloids repel each other. Variants of charge stabilization may include 
the coverage with ionizable molecules which dissociate in water, or the labeling of colloids with charged fluorescent 
markers. For charge-stabilized colloids at large distances, the resulting repulsive force at water interfaces stems from 
a dipole-dipole interaction as shown theoretically for point charges on surfaces [4] and verified experimentally for 
polystyrene (PS) spheres on water-oil interfaces [5]. 

Nonetheless, charged colloids at interfaces also show attractions far beyond the range of van-der-Waals forces. The 
corresponding experimental evidence can be roughly classified as follows, (i) According to Refs. [6-10], PS spheres 
(radii R = 0.25 . . . 2.5 fim) on flat water-air interfaces using highly deionized water exhibit spontaneous formation of 
complicated metastable mesostructures. They are consistent with the presence of an attractive, secondary minimum 
in the intercolloidal potential at distances d/R « 3 . . . 10 with a depth of a few fc^T. The use of water slightly 
contaminated by ions seems to move the minimum further out and to reduce its depth [10, 11]. (ii) In Rcf. [12], 
PMMA spheres with radius R — 0.75 //m were trapped at the interface of water droplets immersed in oil. Here, 
the secondary minimum has been measured at a distance d/R = 7.6 and is reported to be surprisingly steep. The 
tentative explanation of these findings given in Ref. [12] invokes an analogue of long-ranged flotation or capillary 
forces which decay cx 1/d. This interpretation was criticized in Ref. [13] (with which the authors of Ref. [12] agreed 
[14]) and in Ref. [15] which both concluded that possible capillary forces in this system are much shorter ranged, i.e., 
cx d~ 7 , but the authors of these references disagree with respect to the sign of this shorter-ranged force. In yet another 
twist of the story, after completion of our work we encountered the very recent Ref. [16] in which the authors claim 
that long-ranged capillary forces cxl/rf caused by the colloidal charges persist for sub-/xm particles. This conclusion 
is based on measurements of the meniscus shape around single glass spheres with radii 200 . . . 300 /im floating at 
water-oil and water-air interfaces. 

Motivated by the experimental data summarized above and the still incomplete theoretical understanding, here we 
undertake a quite general analysis of capillary interactions between two spherical colloids trapped at fluid interfaces. 
We characterize the system by a general stress field which acts on the interface, e.g., due to a discontinuous electrostatic 
field at the interface, and by forces on the colloid of gravitational and/or electrostatic origin. Special attention will be 
given to the role of a restoring force fixing the interface, e.g., due to gravity or interface pinning. In Sec. II we present 
a free energy model for a single colloid trapped at an interface in the limit of small stresses and forces. The general 
solution for the interface deformation will be used in Sec. Ill in order to determine the effective potential between 
two colloids within the superposition approximation. In view of the differing theoretical results in the literature, the 
derivation of the interface deformation and the resulting effective potential is presented in detail in order to reveal 
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properly the subtleties involved. In Sec. IV we apply the results to the case of charged polymeric spheres on water 
interfaces. It turns out that the constraint of approximate mechanical isolation of the experimental systems renders 
the capillary interaction basically short-ranged. Long-ranged attractive forces oc 1/d only arise through a restoring 
force acting on the interface (which is, however, weak) or in the presence of an external electric field. We shall discuss 
the relation to previous theoretical results, especially in view of the experimental results reported in Rcfs. [12, 16]. 
Directions for further research will be pointed out. In Sec. V we summarize our results. 



II. EQUILIBRIUM STATE OF A SINGLE COLLOID 

In this section we consider as a first step the equilibrium state of a single colloid of radius R at the interface between 
two fluid phases denoted as 1 and 2. The contact angle formed by the interface and the colloid surface is given by 
Young's law, 

coa 0=li^Jl (1) 

7 

where 7 is the surface tension between phases 1 and 2, and ji is the surface tension between the colloid and phase 
i. As a reference state, with respect to which changes in free energy will be measured, we take a planar meniscus 
configuration with the colloid at such a height h lc f that Young's law is satisfied (Fig. 1). This corresponds to the 
equilibrium configuration of an uncharged colloid at the interface if its weight can be neglected — which for generic 
cases is a safe approximation for R < 1 ^im [17] (see also Sec. IV). We model the colloid as a smooth sphere so that 
the system is invariant under rotations around the colloid axis perpendicular to the reference planar meniscus. The 
presence of charges induces a shift of the system (colloid and interface) with respect to the reference state. Here we 
neglect corresponding changes in the surface tensions 7 and 7$; this approximation is expected to be valid provided 
the concentration of charges is sufficiently small. This shift is characterized by the meniscus profile u(r) relative to the 
planar configuration and by the height h of the colloid center. In the reference configuration, the charge distribution 
is assumed to be already in equilibrium. In the following we do not consider the degree of freedom "charge density 
field" explicitly but take it to be fixed to that of the reference configuration. This amounts to neglecting the feedback 
of the interface displacement on the charge distribution. It turns out to be useful to introduce the radius r of the 
three-phase contact line and the angle £ as auxiliary variables (see Fig. 2). 



A. The free energy 

In this subsection we formulate a free energy functional for the degrees of freedom u{r) and h. As shown later, for 
the cases of interest here the deviations from the reference configuration are small enough to justify a perturbative 
treatment, and the free energy can be expanded up to quadratic order in Ah := h — h le f and u(r). We denote with 
II(r) the vertical force per unit area acting on the meniscus surface in the reference configuration. In the case of a 
charged colloid, IT(r) is given by the zz-component of the difference of the Maxwell stress tensor right above and below 
the meniscus, plus the pressure difference acting across the meniscus (incluiding an imbalance in osmotic pressure 
due to the different concentration of ions just above and below the meniscus) [15, 18]. We introduce the following 
dimcnsionless parameter to measure the relative strength of this force: 

e n := f dAIl = —^—[ dr rll(r), (2) 

27T7r ,rcf J5 mc „, r cf 7»*0,ref Jr , Iaf 

where the integral extends over the flat meniscus S'men.ref (the plane z = with a circular hole of radius ro, re f)- We 
assume H(r — > 00) <~ r~" with n > 2, so that the integral converges. In the same spirit we introduce the total vertical 
force F = Fe z acting on the charged colloid in the flat reference configuration, which includes gravity, the electrostatic 
force, and the total (i.e., hydrostatic and osmotic) pressure exerted by the fluids. This leads to the definition 

F 

s F := . (3) 

27r7r , ro f 

These two dimensionless parameters will appear naturally in the course of the calculations. If en and ep vanish, the 
reference configuration is the equilibrium state; it is the global minimum of the free energy functional given in Eqs. (4) 
and (13) below if II = 0, F = 0. The aforementioned perturbative expansion of the free energy can be rephrased as 
an expansion in terms of the small parameters en and ep. 
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FIG. 1: Geometry of the reference state. The equilibrium contact angle fixes the height ft re f 
the contact radius ro, rc f = RsaiO, and the auxiliary angular variable £ re f = 6. 
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FIG. 2: Description of deviations from the reference state. u(r) is the meniscus profile and h is the height of the colloid center. 
The angle £ and the radius ro = i?sin£ of the three-phase contact line are auxiliary variables, which depend on u(r) and h 
through the geometrical relationship h = u(ro) — i?cos£. In this mesoscopic description ro is defined as the position where the 
meniscus profile intersects the surface of the sphere. Ai is the surface area of the colloid exposed to phase i. 



The free energy T of the colloid expressed relative to the aforementioned reference configuration consists of five 
terms: 



•J~ J~ cont ~t~ J~\ men ~t~ vol ~t~ Winter "T" J~\ coll- 

In the following we discuss each contribution: 



(4) 
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• Fluid contact of the colloid. With Ai denoting the surface area of the colloid which is in contact with phase i, 
the surface free energy of the colloid due to its exposure to the phases 1 and 2 is 

•^cont = 71^1 + 72^2 - (7l^-l,ref + 72^2,ref)- (5) 

In Appendix A we express this contribution as a function of Ah and u(r). The final result, valid up to corrections 
of at least third order in e n or £f, reads 

^"cont ^ 7T7N r o,rof) ~ Ah} 2 + irj(r 2 - ro jref ). (6) 

• Change of the meniscus area. The free energy contribution due to variations in the meniscus area relative to 
the reference state reads 

JT men = 7 I dA v/l + |Vw| 2 - 7 j dA, (7) 



where S mcn is the surface of the fluid interface projected onto the plane z = (in which the reference interface is 
located), and V is the gradient operator on the flat reference interface. For small slopes (|Vw| <C 1) one obtains: 

fmcn 7 / dA+ \ 1 dA\Wu\ 2 

>5'iiien\'S'men,ref S'men 

^ Kl{r 2 , cf -rl)+ 1 - 1 f dA\Vu\ 2 . (8) 

Since the w-dependent term is of second order in u, we have approximated the integration domain S men by 
Smen.ref; the corrections are at least of third order in the small parameters en or ef- The first term in this 
expression represents the change in the area of the meniscus which is cut out by the colloid, and in Eq. (4) it 
cancels the second term of Eq. (6). 

• Volume forces on the fluids. We consider the case that the only volume force acting on the fluid phases is 
gravity. The electrostatic forces are active only at the surfaces, where a net charge can accumulate. If the 
spherical colloid is replaced by a long cylinder of radius ro. rc f (Fig. 3) the change in gravitational potential 
energy relative to the reference state due to the displacements of volumes of the fluid phases can be calculated 
easily: 



^voi = ^7 / dA — Y = 7T7 / drr — \L , (9) 



assuming that the mass density of phase 2 is larger than the one of phase 1, i.e., Q2 > £>i; A denotes the capillary 
length 



A = \/i — ^—r (10) 

V \P2 -Pi\9 

in terms of the acceleration g of gravity. For a spherical colloid T vo \ in Eq. (9) has to be corrected to account 
for the specific dependence of the displaced fluid volume on the colloid height h, yielding a slightly cumbersome 
expression [19]. This correction is usually numerically small and we neglect it for reasons of simplicity. As it 
will be shown in Subscc. II C this is justified because the results of interest here are insensitive to the precise 
form of J^oi in the limit A — > oo (A w 1 mm, which is much larger than any other relevant length scale of the 
system). 

• Force on the fluid interface. The aforementioned surface force density II(r) acts on the fluid interface between 
phase 1 and phase 2. The free energy change due to the ensuing displacements of the meniscus is 

Winter --/ dAUu^-f dAnu = -2irf dr rU(r)u(r). (11) 

Here II(r) is the surface force in the reference configuration (II > corresponds to a force pointing upward). 
Changes in the force induced by meniscus deformations and colloidal displacements contribute terms of higher 
orders in e n or £_f in Eq. (11). The replacement of S'mon by Smen.ref m Eq. (11) introduces terms of higher order, 
too. 
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FIG. 3: Long cylinder of radius r 0i rcf immersed vertically into the fluid interface between the phases 1 and 2. 



• Contribution from the colloid. The free energy change due to a vertical displacement of the colloid is 

^coii -FAh, (12) 

where F is the vertical force on the colloid in the reference configuration. Like for Winter, changes of F due to 
deviations from the reference configuration contribute to higher order terms. 

In conclusion, by adding Eqs. (6), (8), (9), (11), and (12) we obtain the following approximate expression for the 
total free energy, which is correct up to second order in e n or e f , and which is a function of Ah and a functional of 
u(r): 



T ~ 27T 7 



poo 




fdu\ 


/ dr r 


\ 




^0, re! 







u 2 1 
2* 



7r 7 [u - Ah] 2 - FAh, 



(13) 



where u = it(r , rof ). 



B. Minimization of the free energy 

The equilibrium configuration of the system minimizes the free energy expression (13). The minimization proceeds 
in two stages. First we seek the minimum with respect to Ah at fixed u(r): 

dT 

= => Ah = u - s F ro,ref, (14) 



d(Ah) 



where we have used the definition (3). Using repeatedly the definitions of h and ro in terms of the auxiliary angle £ 
(Fig. 2) we can compute the change of the contact radius: 

Ar = r - r . ro f = i?(sin£ - sin#) ~ (9 - £)h le{ ~ — — — /i ref = -s F h Te {, (15) 

ro.ref 

plus corrections of second order in sn, £f- In the second step, we minimize Eq. (13) with respect to u(r) at fixed Ah. 
This is a problem of variations with a free boundary condition at r = ro, re f [20]. Variation with respect to u(r ^ ro, re f) 
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yields a second-order ordinary differential equation, 



<Pu Idu 1 1 , x 

TT + ~T" ~ T2 U = 11 r )' 

dr z r dr \ l 7 

while variation with respect to w provides a boundary condition: 

u — Ah 



du 
dr 



»*0,ref 



(16) 



(17) 



where the last equality follows from Eq. (14). The second boundary condition one has to impose on Eq. (16) follows 
from the requirement that the meniscus is asymptotically flat far from the colloid (assuming that U(r — > 00) = 0): 



lim u(r) = 0. 



(18) 



Equation (16) describes mechanical equilibrium of the interface such that the Laplace pressure balances the forces 
acting on the interface. The boundary condition (17) expresses mechanical equilibrium of the colloidal particle: At 
the contact line (i.e., the circle with radius ro at z = uo) the interface exerts a force onto the colloid which has a 
non-vanishing contribution only in z-direction with the magnitude 2njro sin[arctan(ud)] s=s 27T7ro,refwd- This contact 
line force is balanced by the total force F. 

The solution of Eqs. (16-18) can be written in terms of the modified Bessel functions of zeroth order: 



(19) 









(5) 
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where the integration constant A is determined by the boundary condition (17). 



C. Asymptotic behavior in the limit A — ► 00 

For typical values of the parameters, A is of the order of millimeters and therefore much larger than any other 
length scale occurring for experiments with sub-micrometer colloids. In order to study the intermediate asymptotics 
('"o.ref 7 r "C A) of u(r) as given by Eq. (19), we insert the asymptotic expansions of the Bessel functions [21] as A — > 00 
and retain those terms which do not vanish in this limit. Assuming that H(r — > 00) decays sufficiently fast, Eq. (19) 
reduces to 

CA 1 f+°° s 

u(r) ~ r . ro f(en - £f) hi / ds s U(s) In -, (20) 

r 7 Jr r 

where C = 2e~ 7E ~ 1.12 and 7e is Euler's constant. In Eq. (20), we have expressed the integration constant A 
appearing in Eq. (19) in terms of zp by using the boundary condition (17). The first term in Eq. (20) is a solution 
of the homogeneous part of Eq. (16) (with A -1 = in the equation), demonstrating that the limit A — > 00 is singular 
as long as £n 7^ £_f- The second term corresponds to a particular solution of the inhomogeneous differential equation. 
If the surface force II(r) decays algebraically, II (r — > 00) cx r~™, this term decays like r 2 ~ n (since we have assumed 
n > 2), so that the logarithmic contribution is dominant. If n < 2, the asymptotic behavior is no longer given by 
Eq. (20) but a different dependence on A arises. 

At distances r of the order of A, the expression (20) is not valid and a crossover to the exact solution (19) takes place 
in order to satisfy the boundary condition at infinity (Eq. (18)). Figure 4 sketches the behavior of the meniscus profile 
u(r). As can be checked directly in the differential equation (16), one has u(r) ~ (A 2 /7)II(r) cx r~ n asymptotically 
as r —> +00, which expresses the balance between the surface force II(r) on the meniscus and the gravitational force. 
There is, however, another intermediate regime r>A but not too large, in which u(r) decays oc exp(— r/A) and which 
corresponds to the solution of Eq. (16) with II set to zero, i.e., when gravity is balanced by the Laplace pressure 
induced by the meniscus curvature. 

These conclusions, as well as the functional dependence of u(r) on II(r), are robust and independent of the details 
of the boundary condition at r — > 00. In order to support this statement, we consider two cases which implement the 
distant boundary condition differently: 
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FIG. 4: The meniscus solution for the parameter choice A/ro, rc f = 10 2 , ef = 1CT 2 , and a dipole-like stress field ro, rc f H(r)/j = 
0.08 (ro,rcf /r) 6 (en = 2-10~ 2 ). The solid line represents the solution u(r) given by Eq. (19). The dashed line is the intermediate 
asymptotic solution given by Eq. (20). Note that the capillary length A is typically of 0(1 mm). In the present context we 
focus on the length scale ro, rc f, r«A, for which Eq. (20) holds. 



• Pinned interface. In the absence of gravity the interface is assumed to be pinned at a finite distance L from the 
colloid. This corresponds to setting A -1 = in Eq. (16) and replacing Eq. (18) by the new boundary condition 
u{L) = 0. In this case the solution of Eq. (16) is given by 

L 1 [ L s 

u(r) = ro,ref (£n - £f) In ds sU(s) In-, (21) 

r 7 J r r 

where 



:,i := — !— [ drrn{r), (22) 



in analogy to Eq. (2). In the intermediate asymptotic regime (r re f , r -C L), the meniscus profile is then given 

by 



L 1 f + °° s 

u(r) ~ r ,ref (en - £f) In / dssll(s)ln-, 

r 7 Jr r 



(23) 



assuming that IT(r) decays sufficiently fast so that the integrals in Eqs. (21) and (22) converge as L — > oo. 
Equation (23) resembles Eq. (20) with CA replaced by L. 

Pinned curved reference interface. In some experiments the interface between phases 1 and 2 is in fact closed, so 
that the colloidal particle lies at the surface of a large nonvolatile spherical droplet of phase 2 which is immersed 
in phase 1 [12] (Fig. 5) and fixed by certain means (e.g. by a glass plate). The free energy functional (13) 
has to be modified to account for the curvature of the reference interface as well as of the constraint that the 
droplet volume remains unchanged under deformation. To determine the interface deformation, we minimize the 
functional and employ the boundary condition that the droplet is fixed at some point far from the colloid. The 
mathematical details and the corresponding solution for u(r) are presented in Appendix B. Here we quote only 
the intermediate asymptotic behavior (ro,ref ; r "C -Rdrop, where i?drop is the radius of the undeformed droplet): 

u(r) ~ r , rc f (e n - e F ) In CRdl °P _ 1 f dssll(s)\n-, (24) 

r 7 Jr r 

with C a numerical constant given by Eq. (B23). Again Eq. (24) closely resembles Eqs. (20) and (23). 

The physical reason for the occurrence of the singularity in the limit A, L, i?drop - * +oo is that a "restoring 
force" far from the particle is required to yield a well-defined unperturbed interface which allows one to determine 
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FIG. 5: Colloid at the surface of a droplet of phase 2 immersed into phase 1. The radius i?dro P of the droplet, which is spherical 
without the colloid, is usually much larger than the colloid radius R. (The deformation of the droplet has been exaggerated.) 



the deformation u(r) unambiguously. For example, if one takes the limit A — > +oo in the Young-Laplace Eq. (16), 
it is inconsistent to impose the boundary condition u(r — > +oo) = in the corresponding solution. The special case 
= £f, however, is not singular; this will be discussed in Sec. IV. Furthermore, the comparison of Eqs. (20), (23), 
and (24) demonstrates that the functional form of the intermediate asymptotic behavior is independent of how the 
restoring force is implemented. This corresponds to the so-called intermediate asymptotic behavior of the second 
kind [22] , characterized by the following features: 

1. There is a length scale (A, L, or i?drop) which is much larger than the other length scales of the system under 
consideration and which seems — at first sight — to be irrelevant. 

2. Nevertheless, this length scale determines the dominant logarithmic (or more generally, the power-law) depen- 
dence. 

3. The detailed physical origin of this length scale (in the examples wc have considered, gravity, pinning of a 
reference flat or curved interface) does not matter. 

Well known examples of this kind of asymptotic behavior are critical phenomena in phase transitions [23]. In that 
case, it is the microscopic length scale given by the amplitude of the correlation length which cannot be set to zero 
although it is much smaller than the correlation length itself. This microscopic length scale is required to formulate 
the power-law behavior of certain properties of the system, but its detailed physical origin is unimportant for the 
universal decay exponents of the power laws. 

III. EFFECTIVE INTERACTION POTENTIAL OF TWO FLOATING COLLOIDS 

In this section we consider the equilibrium state of two identical colloids floating at the interface at a fixed lateral 
distance d and compute the effective interaction potential V mcn (d) generated by the meniscus. The free energy can be 
derived along the same lines leading to Eq. (13), but with due account for the fact that in the presence of two colloids 
the meniscus slope no longer exhibits rotational symmetry. Nonetheless V mcn depends only on the distance d between 
the centers of the two spheres. The reference configuration is that of two colloids floating on a planar interface with 
the corresponding reference free energy being independent of d. In this case one has for the free energy relative to 
that of the reference configuration a contribution J- mcn + ^vol + Winter from the meniscus and a contribution of the 
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FIG. 6: Top view (projection onto the plane z = 0) of the reference configuration with two colloids, d is the distance between 
the colloid centers. Si and 52 are disks of radius r , ro f, the corresponding circumferences (counterclockwise) are dSi and dS2 ■ 
The projection of the interface is S m cn,rcf = R 2 \(Si (J £2). 



form ^"cont + ^coii from each colloid. (The total free energy includes also the direct interaction between the colloids; 
this contribution will be considered in Subsec. IVB.) From Eqs. (A7), (8), (9), (11), and (12) one obtains: 



T ~ 7 



/ 

JR 2 \(Si 



dA 



2 |VM| + 2^-^ U 



E 

a=l,2 



7 



2ro,ref /as, 



d£ [Ah D 



FAh r 



(25) 



Here, u is the meniscus profile in the presence of two colloids, Ah a are the corresponding heights, II is the vertical 
force per unit area on the meniscus in the reference configuration, and F a is the force on each colloid. By symmetry, 
one has Ah\ — A/12 and F\ = F2. S a are the circular disks delimited by the contact lines of the colloids in the 
reference configuration, dS a are the contact lines, with the convention that we trace them counterclockwise, and 
^mcn.rcf — K 2 \(SiUS 2 ) (see Fig. 6). 



A. Minimization of the free energy within the superposition approximation 

The equilibrium configuration is the minimum of the free energy given by Eq. (25). The minimization procedure 
follows closely Subsec. II B. First, minimizing with respect to Ah a at fixed meniscus height u leads to the height of 
the colloids, 

dT 

~, A , x = => A/i a = u a - e p r , re f , (26) 



where 
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is the mean height of the contact line. In the next step, we minimize with respect to u at fixed Ah a . Variation in the 
interior of the domain M 2 \(S' 1 [J S 2 ) provides a second-order partial differential equation, 



W 2 u- — u= — n, 
\ z 7 

while variation at the boundary dS\ (J 8S2 provides the following transversality conditions [20] : 



du(r) u(r) — Ah a 
dn a r ,rof 



= £p + 



u(r) 



u c 



fO.ref 



r = (x,y) e dS a , 



(28) 



(29) 



where the last equality follows from Eq. (26). In this expression, d/dn a is the derivative in the outward normal 
direction of dS a . (In this way, the triad (e„, e t , e z ) is right-handed, where e„ is the unit vector in the outward normal 
direction, e t is the unit vector in the counterclockwise tangent direction, and e z is the unit vector in the positive 
^-direction.) When applying the transversality condition, in the context of Gauss' theorem one must keep in mind 
that the boundary of the region ]R 2 \(S'i (J S2) consists of the contours dS a traced in clockwise direction with the 
normals directed towards the interior of S a . (The boundary at r — > 00 docs not contribute due to Eq. (30).) In the 
special case of a single colloid, rotational invariance reduces Eq. (29) to Eq. (17). Finally, one has the additional 
boundary condition 



lim u(r) = 0. 



(30) 



The solution of Eq. (28) with the boundary conditions given by Eqs. (29) and (30) is a difficult task. (We are 
only aware of the — already very involved — solution of the homogeneous Eq. (28), i.e., II = 0, with the simplified 
boundary condition du/dn a = e f [24].) For the present purpose, one can use the so-called superposition approximation 
[19, 25, 26], which yields the correct solution in the asymptotic limit of large separation d^> R between the colloids. 
Let u a denote the equilibrium meniscus profile as if colloid a was alone, with II Q and F a denoting the corresponding 
forces. The superposition approximation then reads: 



u 

n 

F 



Ul - 

III 

Fx 



U2, 

-n 2 , 
--f 2 . 



(31) 



Notice that the fields u a (r) and II Q (r) arc defined in the domain R 2 \S Q , while the fields u(r) and n(r) are defined 
in the smaller domain M 2 \(S'i (J S 2 )- Equations (28) and (30) are fulfilled by this approximate solution, but the 
boundary condition (29) is violated: using Eqs. (31) and the boundary condition in Eq. (17) for the single-colloid 
solution, one obtains 



du(r) 
dni 



u(r) — ui 



du 2 (r) 
dni 



»*0,ref 



u 2 {r) 



27rr ,ref JdSt 



d£ U2 



(32) 



and a similar expression for the other colloid with the indices 1 and 2 interchanged. In general, this expression does 
not vanish as required by Eq. (29). If d is large, Eq. (32) can be evaluated by expanding u 2 into a Taylor series around 
the center of colloid 1, yielding to lowest order (ei is the outward directed normal unit vector of dSi) 



du(r) 
dm 



u(r) — ui 



T"0,ref 



redS! 



- jT"o,ref[eiei : VVu 2 (rf) + V u 2 (d)} 



(33) 



in dyadic notation. Inserting the single-colloid solution given in Eq. (20), one finds that this expression decays like 
d~ 2 if ef 7^ En, an d like U(d) ~ d~ n if ef = £n- 

The superposition solution can be used to determine the vertical displacement according to Eq. (26): 



Ah a ^ Ah + u, 

where Ah is the relative vertical displacement of an isolated colloid (Eq. (14)) and 

1 



u : = 



dl U2 



27rr 0i ref jQSx 

is the average of the single-colloid meniscus height at the contact line of the other colloid. 



(34) 



(35) 
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B. Effective interaction potential 



The meniscus-induced effective potential between the two colloids (without their direct interaction) is defined as 

V mcn (d) :=f-T x -T 2 , (36) 

where T is the free energy of the two-colloid equilibrium configuration (Eq. (25)) while T\ = T 2 is the single-colloid 
equilibrium free energy (Eq. (13)). As noted before the energy of the reference configuration is independent of the 
separation d and drops out from Eq. (36). We insert Eqs. (31), (34), and (35) into Eq. (36) and exploit the invariance 
of the free energy under exchange of the colloids, i.e., the symmetry under exchanging indices 1 <-> 2. After carrying 
out some analytic manipulations one finds the following expression for the effective potential: 



V men (d) 



dA 

\(Si[JS 2 ) 



J 

Ssi 

+ — — I di[u- u 2 f - 2Fu . 

nVef JdSi 



- / dA 

'Si 



7 

7(Vwi) • (Vu 2 ) + -yU^ - 2II1U2 



(37) 



The first integral accounts for the change in surface energy, gravitational potential energy, and surface-stress potential 
energy of the meniscus due to the overlap of the meniscus deformations caused by the two colloids. The second integral 
is the corresponding change due to the fact that the interface is reduced by an amount S\ compared to the single- 
colloid case because of the presence of the second colloid. The third integral is the change in surface free energy of one 
colloid due to the extra meniscus deformation induced by the second colloid. The last term is the change in energy 
due to the vertical displacement of one colloid by this extra meniscus deformation. 

For the mathematical manipulations to follow, it is suitable to rewrite Eq. (37) by applying Gauss' theorem to the 
integrals involving V« and by using the fact that the functions u a fulfill Eq. (16) individually, e.g., 



7 / dA (Vtti) • (Vu 2 ) 

Jr 2 \(SiUSh) 



dA [V • (u 2 Vui) - u 2 V 2 ui] = 
\(Si U s 2 ) 

d(uiu 2 ) f 



- 7 f dt 

IdSi 



dni 



/ dA 

Ir 2 \(S! u s 2 ) 



II1M2 - -^UlU 2 



(38) 



Thus one obtains from Eq. (37) 



dATl 2 u 2 

JSi 



V men (d) ~ - / dATl lU2 + 

JS. 2 \(Si\J S 2 ) 

JdSi om 2 J dSl dnx 
+ I dl\u- u 2 f - 2Fu. 

ro,ref JdSi 



(39) 



In this form all the integrals, except the first one, are performed over bounded domains. This allows one to carry out 
a Taylor expansion which yields a uniformly valid asymptotic expansions of the terms as d — > 00. In the following 
calculations, the integrals are evaluated in polar coordinates with the origin at the center of S a and r a is the position 
vector with respect to this origin. In particular, r 2 = de x is the center of colloid 1 with respect to the center of colloid 
2, and 17 = — de x is the center of colloid 2 with respect to the center of colloid 1 (Fig. 6) so that, e.g., u a = u(|r Q |). 



1. In order to compute the leading asymptotic behavior of the first term in Eq. (39) as d - 
cases: 



00 we distinguish two 



(a) £n 7^ £f- In this case, u ~ logr (Eq. (20)) while II 
contribution stems from the region near Si and 



r " (n > 2) as r — > 00. Asymptotically the main 



dA ITi u 2 ~ u(d) 
! \(SiUs 2 ) 



dA IIi = 27T7r ,ref en u(d), 

! \Si 



(40) 
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employing the definition (2). It will turn out that there is no need to compute also the next -to-leading 
term in order to obtain V mcn (d — > oo) as the leading term will be part of the leading contribution to V men 
and it will not be cancelled by other contributions. 



(b) £n = £f- In this 



case u <~ r 



if n(r 



The regions contributing mainly to the integral 



are the neighborhoods of Si and S 2 . We employ the Taylor expansion of the integrand (which is valid up 
to a maximum order depending on how fast II decays) to evaluate the leading and the next-to-leading 
contributions if n > 4 (using dyadic notation): 



u 2 + ri ■ (Vtt)j + 2 VlVl : ( VVu )2 + • • ■ 



/ dAIliu 2 ~ / dATl 

Jr 2 \{Si\J S 2 ) JR 2 \Si 

+ [ dAu 2 [Ui + ...} ri= _ dex 

Jm 2 \s 2 

1 f°° 
~ 27r7r ,ref £n + -7rV 2 u(<i) / dr r 3 U(r) 

J rn rof 



r 2 =de x 



+ n(d) 



rfA M 2 + • • • • 

! \s 2 



(41) 



We can simplify the result further by evaluating the last integral by repeated partial integration with the 
explicit solution for u given in Eq. (20): 



1 



dA u(r 2 ) = -7rr f 

R 2 \5 2 2 



£n 



2u 



r"0,ref 7 r 0,ref Jr 0jrof 



dr r 3 II(r) 



One finally obtains 



/ dA 111 u 2 ~ 27r7r . ro f £n 

Jr 2 \(SiU s 2 ) 



(42) 



(43) 



+ 



u(r 2 ) 

R 2 \S 2 



1 



2mq 



11(d) + 



where we have employed Eq. (28) (with A 1 = for simplicity). Convergence of the integral of u imposes 
the more stringent constraint n > 4 on the decay of II (see Eq. (42)). 

The validity of formulae (40) and (43) has been checked explicitly by comparing their predictions with the 
numerical evaluation of the integral for a surface stress of the form II oc r~ n with n > 4. 

2. The second term in Eq. (39) can be estimated in the limit d — > oo by expanding the integrand into a Taylor 
series, which is uniformly valid in the integration domain: 



/ dA n 2 u 2 ~ / dA [II 2 u 2 + ... } r2=de ~ Trrg rcf 11(d) u(d) + ... 

JS! JS! 



3. The third term in Eq. (39) reads: 

d(uiu 2 ) 



dl . 
as! dn i 



271 



dip [n • V(uiu 2 )] 



(44) 



(45) 



* iy { 



(n • V) 



Mi ( u 2 + ri • (V«) 2 + 



-riri : (VVu) 2 



~ 27rr , ro fu((i) 



du 
dr 



r 2 =de XJ , ri=ro ref 

^r ^ 2 u{d) 



2irr a , Ici e F u(d) - ^ r o,r C f [ £ f 



2mq 

>"0,ref 



dr r=r 0jrcf 

n(d) + . . . 



(46) 



using the boundary condition for u in Eq. (17). 
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4. The fourth term in Eq. (39) is evaluated analogously: 

,2 /-27T 



I dlpl = f\ [ ri . Vu 2 ] (47) 

J 0Sl dm 7o L 2Jri=ro -' v 7 

* j\ {r 1 -[(Vn 2 ) 2 + r 1 .(VV W 2 ) 2 + ...] r2=de J j 



1 ri=r ,ref 

„2 V72..2; 



~ 7rr$ iref VV(d) + ... . (48) 
5. The fifth term in Eq. (39) is given by 

1 f f 2rr 

dl [u2 — u} 2 = / dtp u 2 — 2nu 2 



fO.ref JdS! 



^rl rcf [V 2 u 2 (d)-2uV 2 u(d)]+... 



2 

1 + . . . . (49) 



1 2 

^O.ref 



V 2 u 2 (d) + ?-n(d)u(d) 



6. In order to evaluate the sixth term in Eq. (39), we expand U2 in the definition given by Eq. (35): 

u = I dlu 2 ~ u(d) + -r 2 , cf V 2 u{d) + ■ ■ ■ ~ u(d) - -^r 2 rof 14(d) + . . . . (50) 

27rr , re f J dSl 4 ■ 47 

The asymptotic behavior of the effective potential is finally obtained from Eq. (39) by collecting all terms. There are 
two qualitatively different cases: 



£n 7^ £f- The limit A — > 00 is singular and the single-colloid meniscus profile exhibits a logarithmic dependence. 
One finds that the leading contribution is provided by the terms proportional to u(d): 

CX 

V mcn (d) ~ 27T7r ,ref (£f - £n) u(d) ~ - 27r7ro ;I . ef (ef - £n) 2 In — , {R < d < A), (51) 

which represents a long-ranged attractive effective potential, irrespective of the sign of the forces F and II acting 
on the system. 

gn = £f - The single-colloid meniscus profile decays like u(d) ~ <i 2_ ™ if 11(d) ~ c? _ ™. The leading contribution 
is proportional to 14(d), because the terms proportional to u(d) cancel each other: 



V men (d) ~ - 2 n(d) / dA u, (i? < d < A). (52) 

JR 2 \S 

This correspond to a shorter-ranged effective interaction, which in principle can be either attractive or repulsive 
depending on the precise form of the function II(r). In the particular case that II(r) decays monotonically to 
zero, e.g., II(r) cx r _ ", it is easy to check that V mcn amounts to a repulsive force, because the sign of u is 
opposite to the sign of 14. 

We have seen that the error committed by the superposition approximation in satisfying one of the boundary 
conditions, Eq. (33), decays like 14(d), too. This suggests that the corrections to the superposition approximation 
could modify the precise value of the constant factor acting as an amplitude in Eq. (52), and a priori it cannot 
be excluded that there are cancellations leading to a vanishing amplitude, and therefore to an even faster decay 
for large d. Thus the superposition approximation might not be reliable enough for calculating V mcn if en = £f- 

If one traces back the origin of the dominant contributions to Vin 0n , one finds that in both cases only the first, 
third, and sixth term in Eq. (39) are relevant. They correspond physically to the effect of the overlap of the two 
single-colloid meniscus profiles and the change of the colloid height (first integral and the term —2Fu in Eq. (37), 
respectively). 
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IV. APPLICATIONS AND DISCUSSION 



Equations (51) and (52) describe the asymptotic behavior of the meniscus-induced effective intercolloidal potential 
and thus represent a central result of our analysis. They provide the explicit functional dependence on an arbitrary 
stress field II(r) which decays sufficiently fast. The assumptions entering their derivation are (i) that the deviations 
of the meniscus profile from the reference configuration arc small, allowing one to confine the analysis to a free 
energy expression which is quadratic in the deviations, and (ii) the superposition approximation, which expresses 
the two-colloid equilibrium state in terms of the single-colloid state. The analysis shows that the limit A — > +oo is 
non singular only in the case ep — £n- 



A. Flotation force 



Equation (51) can be used to determine the flotation force. There are no stresses acting at the meniscus, II = 0, 
while the force F on the colloids is due to their weight corrected by the buoyancy force. Accordingly, Eq. (51) reduces 
to the flotation potential 

V Rot (d) = 27T7r ,ref£F u(d) (53) 

= 2^ 7 Q 2 ln-4 , A), 
o A 

where Q := £_Fro, rc f is known as the capillary charge [17], by analogy with two-dimensional electrostatics. The asymp- 
totic form for d <C A originates from a potential proportional to the modified Bcsscl function Ko(d/\) (see Eq. (19)) 
[25, 26]. The order of magnitude of the capillary charge is Q <~ ro, re f (-R/A) 2 . For a typical value 7 = 0.05 N m _1 at 
room temperature and for colloids with a mass density of the order of 1 g cm -3 , the prefactor of the logarithm can 
be estimated as 

Therefore, compared with the thermal energy, the flotation force is negligible for submicrometer-sized colloids. 



B. Electrically charged colloids 



Another application is the case of electrically charged colloids. If one of the liquid phases is water, the charge of the 
colloid is screened (the Debye length of pure water is w 1 /an and smaller in the presence of an electrolyte), and the 
effective electric field is that of a dipole oriented perpendicularly to the fluid interface [4, 5, 27-29]. The electrostatic 
field decays as r -3 and the stress on the meniscus as II(r — ► 00) oc r~ 6 . Both the electrostatic stress and the osmotic 
ionic pressure decay in the same manner [15]. Thus the total intercolloidal potential at intermediate distances is 

Vtot = a^- + V mcn (a > 0) . (55) 

If gravity is neglected both as a force on the colloid and as a restoring force for the interface (A — ► 00), then one can 
indeed show that the ensuing condition of mechanical isolation (no net force on the system) leads to en = £f, i.e., 
precisely the situation for which the limit A — > 00 is non-singular. To see this, we consider the total stress tensor 
T which consists of the Maxwell stress tensor (due to the electrostatic field) and a diagonal osmotic pressure tensor 
(due to the electrolytes) [18, 30]. At interfaces T can be discontinuous. The total volume V of the system is divided 
into volumes V\, V2, and V3 (see Fig. 7 for the explanation of the notation in the following equation) . The total force 
reads (the superscript +( _ ) denotes evaluation on the positive (negative) side of the oriented surface, i.e., the side the 
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arrows in Fig. 7 point to (do not point to)) 

dA-T = [ dV(V-T)+[ dA-(T+-T-) 

JVi\JV 2 \JV 3 Js men, ref U Si u s 2 

= [ dV(V-T)+[ dA-(T + -T") + 

JVi u v 2 JS men , ref 

[dV(V-T)+[ dA ■ (T + — T~) 

^ V3 JSi (J 5 2 



cL4IIe z + Fe z . (56) 

$men,ref 

In the first line, we have applied Gauss' theorem with due account of the possible discontinuities of the tensor T 
accross the interfaces. In the second line, V • T = in the fluid phases V\ and V2, because the counterion distribution 
of the present reference configuration is the equilibrium distribution and thus locally force free. This distribution is 
considered to be fixed. The second term in the second line is the total force on the meniscus (which can have only a 
normal component), while the terms in square brackets sum up to the force F acting on the colloid. 
Thus the vertical component of the total force is 

dA ■ T = 27r 7 r , rof (e n - e F ). (57) 

ot 

If it vanishes, as it is the case for an isolated system, then en — £f- According to Eq. (51) this implies that the 
long-ranged logarithmic contribution to V m cn is absent and thus the limit A — > 00 is regular. Physically this means 
that there is no need for a restoring force acting on the fluid-fluid interface when the deformation is created by 
localized internal stresses. Instead, according to Eq. (52), one obtains a potential y mon oc d~ e for the present case of a 
dipolar electric field (see above). This shorter-ranged potential cannot counterbalance the direct electrostatic dipolar 
repulsion ~ d~ 3 . Such a counterbalance would be needed for a straightforward explanation of the aforementioned 
experimentally observed attractions. 

The line of argument to explain the absence of the logarithm contribution to V men was already put forward in 
Refs. [13] and [15], where exclusively the case II(r) cx r~ 6 has been addressed. Our detailed analysis complements and 
generalizes these contributions. In Ref. [13], V mC n(d) is estimated by taking into account only the degree of freedom 
"meniscus profile", w(r), and considering only the change in meniscus area due to the superposition of the dimples. 
This corresponds to retaining only the term <x (Vui) • (VU2) in expression (37). Although the type of d-dependence 
obtained that way is correct, the sign of the force turns out to be wrong (attraction instead of repulsion). The reason 
is that the contributions cx IT 7/2 and oc Fu in Eq. (37) are equally important as the retained term. This is taken into 
account in the more detailed analysis of Ref. [15], where the limiting case R = (point-particle) is considered from 
the outset, so that "height of the colloid" is not an independent degree of freedom, i.e., h = u(0). Correspondingly, 
F is set to zero (sf — 0), and its effect is modelled by a Dirac-delta contribution to the stress II(r) such that en = 0. 
The potential V mcn (d) calculated this way corresponds to keeping the three terms mentioned above which are relevant 
(the first and third term in the first integral of Eq. (37) and the last term of Eq. (37)) and to setting A — > 00, since 
this limit is regular when ep = en = 0. Our analysis has shown that the terms which are dropped in the limit 
R — > (second and third integral in Eq. (37), second, fourth, and fifth terms in Eq. (39)) yield indeed a subdominant 
contribution to Eq. (52). However, the integral appearing as a prefactor in Eq. (52) is divergent for R — > 0, so that in 
Ref. [15] a short-distance cutoff a has been introduced which is expected to be of the order of R. The precise value 
of a depends on the details of the implementation of this cutoff; for the example used in Ref. [15], the application of 
Eq. (52) yields a = r 0iro f • 

In the presence of gravity mechanical isolation is violated and as we have shown en — sf oc (R/X) 2 . (The force- 
balance argument (Eq. (56)) can be easily generalized to include the effect of the gravitational volume force.) In the 
case of a curved reference interface (corresponding to the experiment reported in Ref. [12] using colloids trapped at 
the interface of water droplets in oil), Eq. (56) is replaced by 

dA-T = I dAHe r + Fe z , (58) 

t S> men, ref 

where e r is the radial unit vector of the unperturbed spherical droplet. Force balance in the vertical direction yields a 
factor e r e z — cos-0 (see Fig. 5) in the integral, which can be expanded for small angles in the limit i?drop — » 00. This 
leads to a curvature-induced correction en — £f oc (ro. ro f/i?drop) 2 even for mechanical isolation. Thus we see that 
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phase 1 




FIG. 7: (color online) In the reference configuration the whole system is divided into volumes Vi, V2, and V3. Volume Vi 
(enclosed by the upper dashed curve) includes phase 1, volume V% (enclosed by the lower dashed curve) includes phase 2, 
and volume V3 includes the interior of the colloid. The arrows indicate the direction in which the surfaces (including the 
infinitesimally displaced ones) are oriented: Stot encloses the whole system, S , mcn ,rof is the interface between phase 1 and phase 
2, and Sim denotes the interface between the colloid and phase 1(2). 



independently of the details of the implementation of the distant boundary conditions the logarithmic term in the 
potential Kn C n (Eq. (51)) has a strength which is proportional to (colloid size/system length) 4 , which is nevertheless 
too weak to explain the reported attractive total interaction. 

1. The experiment reported in Ref. [12] 

Mechanical isolation is violated in the presence of an external homogeneous electric field Ee z . For the experimental 
setup used in Ref. [12] it cannot be ruled out that such an external field may have distorted the measurements [31]. 
Since this is the only experiment which provides quantitative information about the secondary minimum in Vt t, we 
discuss the case of an external field in more detail. In this experiment position-recording measurements were performed 
on a hexagonal configuration of seven trapped colloids on a water droplet immersed in oil (i?drop ~ 32i? = 24/im). 
The latter was confined between two glass plates and the droplet stuck on the upper glass plate with a contact angle 
close to 7r. Residual charges might have resided on the upper plate [31]. The measurements yielded the position of 
the secondary minimum, d mul 7.6i?, and the curvature at the minimum, Vj." t (c? m i n ) ss 12.94 fc^T/i? 2 . With regard 
to the latter value we remark that systematic corrections can be estimated to lower V^ t by a factor 2 to 3. These 
systematic corrections include (i) multiparticle effects and (ii) center-of-mass movement of the droplet. We have 
estimated the effect of (i) by carrying out Langevin simulations of the seven particle system using the intercolloidal 
potential Eq. (59) below. As for the effect of (ii), any shape deformation induced by the moving colloids changes the 
center-of-mass position of the droplet since the droplet is fixed to the upper glass plate. The corresponding change 
in the gravitational energy of the droplet translates into a weak confining potential for the colloids which limits the 
stochastic movement of the center-of-mass of the seven colloids. This effect might be part of an explanation for the 
absence of center-of-mass movement observed in Ref. [12]. According to Eqs. (51) and (55), the total intercolloidal 
potential in an external field is 

V tot = a^f - bln^- (a,b>0). (59) 
ar a 
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Using the aforementioned experimental data for c? m ; n and V^ t (d m i n ), one obtains from Eq. (59) b w 249 fc^T, a w 
83d^ in and Vtot(^min) ~ —275 fe^T which is surprisingly deep, even if reduced by a factor 2 to 3 due to the systematic 
corrections mentioned above. Furthermore, from Eq. (51) we deduce b — 2n-yrl rcf (en — £_f) 2 and with 7 w 0.05 N/m 
we find |en — £f\ ~ (2 nm/ro jTO f). The long-ranged meniscus deformation (see Eq. (20)) is on the scale of nanometers. 
The short-ranged meniscus deformation near the colloid can only be evaluated with a specific microscopic model for 
Il(r). For a rough estimate of en, we consider the colloid charge to be concentrated on the surface. The asymptotic 
behavior of the stress tensor in this case is given by H(r) = a/(4n) (fcsT)/r 6 [4, 15]. If we assume this form of the 
stress tensor to hold for all r, we find en ~ 10~ 4 / sin 5 9 for the values of 7 and a given above. For not too small 
contact angles, this a posteriori justifies the pcrturbativc approach which we have adopted. 

If the system has a net charge q, then \qE\ = 27r7r 0!TO f |en — £f\ = V^irjb (see Eqs. (57) and (51)). Using the values 
for b and 7 as given above we find 

\qE\ w 3.6 x lO^Vm" 1 . (60) 

For the value q s=a 2 x 10 5 e quoted in Ref. [12], this yields a relatively small electric field E w 1.8 x 10 4 Vm _1 , 
indicating how sensitive the system can be to spurious external fields. Alternatively, an electric field E ~ 10 3 Vm _1 
is sufficient for the meniscus-induced logarithmic potential to have a depth of the order 6 ~ 1 ksT. Thus the external 
field offers the possibility to tune easily the capillary long-ranged attraction and to manipulate the structures formed 
by the colloids at the interface. 

2. The experiment and analysis reported in Ref. [16] 

In Ref. [16], experimental results for the meniscus deformation around glass particles of radii 200. . . 300 /im trapped 
at water-air or water-oil interfaces are reported. The data for the meniscus slope u'(tq) at the contact circle imply 
(Eq. (17)) Ef ~ 0.2 (water-oil interface), about fifteen times larger than the corresponding ef, s caused by gravity 
alone. Furthermore, the reported meniscus shape for one sample contains a logarithmic part which is consistent with 
£f — £n ~ 0.1 3> e_F, g (see Eq. (20)). As in the experiment analyzed in the previous subsection, the experimental 
observations could be understood within the framework of the theoretical model we have developed in terms of an 
external electric field violating mechanical isolation. Yet independently of us, the authors of Ref. [16] have developed 
a theoretical model based on the same physical hypotheses as our approach which, they claim, explains the observed 
long-ranged meniscus deformation. Here we would like to note two important errors which flaw their analysis: 

1. Eq. (3.15) in Ref. [16], in terms of which they interpret their observations, can be obtained by inserting the 
large-distance (r » ro,ref) asymptotic behavior of II(r) in our Eq. (20) (equivalent to their Eq. (3.14) up to 
an additive constant) both into the integral term and into the definition of en (Eq. (2)). Since the dominant 
contribution to the integral defining en stems from points r ~ ro.refj this procedure is clearly inadmissible. As a 
consequence, they obtain a wrong, nonvanishing prefactor of the logarithm, in spite of their explicit consideration 
of mechanical isolation (see their Eq. (6.6)). 

2. In order to calculate the intercolloidal effective potential within the superposition approximation, the formula 
relating V men (d) to the meniscus deformation u(d) in the presence of gravity alone (i.e., Eq. (53)) is applied even 
though 11^0. Thus, an additional term contributing to the leading logarithm is not included (compare their 
Eq. (3.16) and our Eq. (51)). 

In Ref. [16] also a numerical analysis is carried out. A detailed study of the relation between the results of this 
numerical analysis and our theoretical predictions will be published elsewhere [32]. 

It may be possible that the presence of an external field is consistent with the data from Refs. [12, 16], as well as 
with the presence of the secondary potential minimum observed in the experiments using planar troughs, in particular 
in the cases d m ; n > 10i?, which fall into the intermediate asymptotic regime considered here. But this still remains 
as an open problem. 

C. Outlook 

Given that already nanometer distortions of the meniscus produce noticeable attractions, the surface topography 
of colloids might be relevant. In Ref. [10], colloidal surface roughness is proposed as an explanation of the attraction. 
The meniscus contact line is assumed to be pinned at defects on the colloid surface caused by surface roughness. This 
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imposes a different boundary condition for the meniscus at contact, which is then deformed even in the absence of 
electrostatic forces. The corresponding analysis in Ref. [10] is concerned only with the term oc (V«i) • (VM2) in the 
expression (37), leading to the conclusion that V men (d) decays as d and corresponds to an attractive potential with 
a strength of 10 4 ksT for meniscus deformations of the order of 50 nm. This conclusion, however, cannot be simply 
carried over to the case of charged colloids. As we have shown the contributions of the terms oc H1U2 and oc Fu in 
Eq. (37) are relevant and can change the qualitative behavior of V mcn (d) even in the limit of point-colloids. It would 
be worthwhile to generalize the analysis of Ref. [10] along the line of arguments presented here in order to assess the 
importance of surface roughness. This should be complemented by more precise experimental information about the 
actual colloidal topography. 

Recently, an explanation based on a contaminated interface has been advocated [33]. The air-water interface 
would be actually a two-dimensional emulsion consisting of hydrophylic (water) patches and hydrophobic (silicon oil) 
patches. The colloidal particles (hydrophobic in character according to Ref. [33]) would cluster in the hydrophobic 
patches. Thus, confinement of the colloids by finite-size hydrophobic patches would give the impression of an effective 
intercolloidal attraction. At present, this explanation is only of qualitative nature. 

Thermal fluctuations of the interface position around its mean value u(r) also induce an effective interaction 
between the colloids which confine these fluctuations (Casimir-like force) . Using a Gaussian model of the fluctuating 
interface in analogy to the procedure employed in Ref. [34] for calculating fluctuation-induced forces between rods in 
a membrane, one finds for uncharged colloids a fluctuation-induced potential Vfl uc ~ — fcgT (r , ro f/rf) 4 , which is too 
small and falls off more rapidly than the intercolloidal dipole repulsion. Here, the generalization to the charged case 
might give hints for the effective attractions between very small particles trapped at interfaces. Concerning particle 
sizes well below the Debye length, one should modify our analysis to account for the overlap of the screening ionic 
clouds (this would affect, e.g., the superposition approximation for the stress field II, Eq. (31)). 

Finally, in Ref. [35] the attraction of particles trapped at a nematic-air interface is reported and an explanation in 
terms of a logarithmic meniscus deformation is proposed which parallels the explanation given in Ref. [12]: in this case, 
the deformation would be caused by the elastic distortion induced by the particles on the nematic phase [35]. Our 
detailed theoretical treatment shows that no long-range (logarithmic) meniscus distortion can arise on an interface 
if the system is mechanically isolated and the excess free energy of the perturbed interface is correctly described by 
an expression like Eq. (13). Thus it appears that the simple explanation of the observed colloidal attractions given 
in Ref. [35] is not correct. However, it is not clear whether the free energy of a distorted nematic-air interface is 
equivalent to that of a simple fluid interface due to the long-ranged interactions in the nematic bulk caused by defects. 
Thus, a generalization of our theory to interfaces involving nematic phases would be desirable in order to assess the 
possibility of long-ranged colloidal attractions in more detail. 

V. SUMMARY 

We have analyzed the effective force induced by capillary deformation between two smooth spherical colloids floating 
at a fluid interface. The relevant degrees of freedom are the meniscus deformation and the height of the colloids 
(Fig. 2), whose equilibrium values are given by minimization of a free energy functional. This functional was derived 
assuming small deviations from a reference configuration (Fig. 1). It incorporates the surface tensions of the three 
interfaces involved ("colloid - fluid phase 1", "colloid - fluid phase 2", and "fluid phase 1 - fluid phase 2"), the 
potential energy of the colloids under the action of a force F, the potential energy of the fluid interface in an arbitrary 
surface stress field n(r), and the potential energy due to a restoring force acting on the interface (Eqs. (4)-(13)). The 
effective intercolloidal potential (Eq. (36)) was calculated in the limit of large separations by using the superposition 
approximation (Eq. (31)). We have shown in this limit that the contribution to the effective potential by the interfaces 
"colloid - fluid phases" is subdominant. If the total force acting on the system, consisting of the two colloids and 
the meniscus, does not vanish, the presence of the restoring force acting on the fluid interface is essential — although 
its precise form does not matter (Subsec. II C). In this case, the effective interaction is long-ranged and attractive 
(Eq. (51)). If the total force vanishes, the restoring force is irrelevant altogether, the effective interaction is shorter- 
ranged (Eq. (52)), and it cannot be computed reliably within the superposition approximation. 

As an application, we have considered the case of like-charged, micrometer-sized particles when the capillary 
deformation is due to the ensuing electrostatic field. We have discussed how one can tune the long-ranged attractive 
interaction by an external electric field, but we conclude that the experimentally observed attraction in an isolated 
system cannot be explained within the present model. Possible directions for future research such as colloidal surface 
roughness and fluctuations of the interface have been discussed. 
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APPENDIX A: FREE ENERGY OF COLLOID-FLUID CONTACT 



In this appendix we determine the contribution to the free energy due to the exposure of the colloid to the two fluid 
phases. Here we consider the general case of no rotational invariance, so that the meniscus height u(r, ip) depends on 
the distance r from the axis through the colloid center and on the angle of revolution ip around this axis (Fig. 8). 
Accordingly, the auxiliary variables r (ip) and £((/?) also depend on the revolution angle and one has 

r (¥>) = -Rsin£(p) (Al) 

and 



h = u(ro(<p), ip) - Rcos£(<p), (A2) 

where u(ro(ip), ip) is the meniscus height at contact. The surface areas of the colloid in contact with phase 1 and 
phase 2 are 

/■27T r£M /-27T 

A 1= R 2 dip di{>smiP = R 2 dtp[l- cos £(ip)} (A3) 
Jo Jo Jo 

and 

A 2 =4ttR 2 - A u (A4) 

respectively. The expression for the free energy in Eq. (5) is based on these formulae for the special case £ ie i(<p) = 9 
(Fig. 1) and on Young's law (Eq. (1)): 

r 2TT 



•7"cont = lR 2 / dip [COS 2 9 — COS # COS £((/?)] = 

Jo 

= 7,1 r2 I dtp [cos £(<p) - cos 6} 2 + -jR 2 / d<p[cos 2 9- cos 2 £(ip)} 
2 Jo 2 Jo 

1 1' 2t: 1 l' 211 

= - 7 dip [u{r (ip),ip) - Ah} 2 + - 7 dip [rl(p) - rg jref ]. (A5) 



The second term, which arises upon completing the square, represents the change of the meniscus area which is cut 
out by the colloid. Since u and Ah are already of first order in £n or ep, one can replace r (ip) by ro, re f m the first 
term and obtains 



T cont ~ -j J dp [u(r ,rcf , p) - Ah} 2 + -7 J dip [rg (ip) - r^ ref ] , (A6) 

plus corrections of at least third order in en or sp. In the special case of rotational invariance, this expression reduces 
to Eq. (6). For the purpose of Sec. Ill one can rewrite Eq. (A6) as the line integral 

T cont ~ I d£ {[u- Ah} 2 + [rl r 2 ^} } , (A7) 

^0,rcf JdS 

where S is a circular disk of radius ro, re f , here centered at the origin, and dS denotes its circumference. 



APPENDIX B: DEFORMATION OF A SPHERICAL REFERENCE INTERFACE 



In this appendix we calculate the deformation of the interface of a spherical droplet due to a floating colloidal 
particle. The reference configuration corresponds to a colloid floating at the surface of the spherical droplet such that 
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FIG. 8: Geometric description of the three-phase contact line between the fluid-fluid interface and the surface of the colloid. 
Side view (a): The auxiliary variables ro(<p) and in general depend on the angle of revolution <p. ip denotes the polar 
angle of points on the colloid spherical surface. Top view (b): The dashed line is the circumference of radius R of the spherical 
colloid. The solid line is the projection of a noncircular contact line. 



Eq. (1) is fulfilled. Here we will consider only axially symmetric configurations (see Fig. 5). In terms of the polar 
angle ip £ [0, 7r], the distance r measured along the unperturbed spherical droplet surface is r = PdropV'- (Note that in 
this appendix we consider u(-) and II(-) as functions of ip instead of r.) We define ipo :— ro,ref/-Rdrop and we introduce 
the angle ipi with 7r > ipi 3> "00 a t which another boundary condition (to be specified below) holds reflecting the fact 
that the droplet is fixed in space. With this additional boundary condition we model closely the actual experimental 
setup of Ref. [12], where the droplet was attached to a glass plate with a contact angle less than ir [31]. 
For the present case the free energy given by Eq. (13) must be modified in several ways: 

• In spherical coordinates, Tm Cn reads 



•^mcn = 7 J dA (-Rdrop + u) 2 J 1 + 



-Rdrop 



-7 



dA i?drop ' 



(Bl) 



where dA — 2n simp dtp is the differential of the solid angle, and V = e^(d/dip) is the gradient operator on the 
unit sphere. As before, this functional is now expanded for small deviations from the reference configuration, 
and we obtain 



^mcn ^ 7 7r ( r 0,rof - r o) + 27r 7 / # sin V 7 

J 1p 



1 fduV 2 



2-/? c j r0 p'u 



(B2) 



where we have neglected terms evaluated at ip which vanish in the limit i?dro P - * 00 • 



• The term Winter in Eq. (11) has to be amended by the extra work upon deformation which is caused by the 
pressure difference in the reference configuration between the interior and the exterior of the droplet, AP = 
27/i?dro P : 



f dA R 2 drop Hu+ (AP) j du (i? drop + uf 



f4>i 

~ — 27T7 / dtp simp 

J j/> 



d2 

J Mrop 

7 



Ilu+(2R dlop u + 2u 2 ) 



(B3) 



Since AP is of order e°, we have to keep also terms oc u 2 in this expression for consistency. For simplicity we 
have neglected the contribution due to the change of S men during the deformation. Such a term would affect 
the behavior of the meniscus only near ipo and vanishes in the limit Pdrop — * oo. 
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• The change in volume of the perturbed droplet reads 



AV 



dA -(i?drop 

Omen 



u) 3 



-^d r0 p 

<S'men ) ref 



(B4) 



where S is a contribution due to the spherical shape of the colloid; it depends on Ah and Aro but can be safely 
neglected. (This is the same kind of term encountered in calculating !F vo i, see Eq. (9).) Expanding AV for small 
deformations u up to linear order yields 



AV ~ 2nR^ Top / dip simp u. 
J 4> 



(B5) 



Physically, the deformation of the interface occurs under the constraint that the droplet volume remains un- 
changed. This constraint (AV = 0) will be implemented in the free energy functional by means of a Lagrange 
multiplier /i which itself turns out to be linear in the parameters en and Sf- This justifies keeping only the 
linear term in Eq. (B5). 

The contribution .Fcont (see Eq. (6)) is concerned only with quantities at tpo, the colloid surface, and thus in the limit 
of large droplet radii this contribution are the same as in the planar case. The free energy related to the work done 
upon moving the colloid, T co \\ (see Eq. (12)), remains also unchanged. For simplicity we also neglect gravity and set 
A -1 = 0. In conclusion, the free energy functional to be minimized is the sum T conil + ^mm + Winter + -^coii phis the 
constraint term (27T7/i?dro P ) M (where fi is a dimensionless Lagrange multiplier): 



T = 27T7 / dtp sin ip 

ho 



du 



R 2 

liRdio P u -^IIu 

7 



7r 7 [u - Ah} 2 - FAh. 



(B6) 



The quantity — 7/i/i?drop can be interpreted as a homogeneous pressure field enforcing the constant -volume constraint. 
We note that the terms linear in u present in Eqs. (B2) and (B3) have cancelled in the total free energy. Minimization 
with respect to Ah yields the same result as in Eq. (14). Subsequent minimization with respect to u(ip) leads to 



sin-0^-) + 2u = ^n + ^drop 

sin ip dip \ dip J 7 



and (compare with Eq. (17)) 



du 
~ckp 



— £fR& 



rop ■ 



The solution must satisfy the constant-volume constraint 



dip sinip u(ip) = 0. 



(B7) 



(B8) 



(B9) 



ipo 



Finally, the second boundary condition mentioned at the beginning expresses the physical requirement that the droplet 
is fixed in space. (Otherwise, the application of a localized force at the interface would shift the droplet as a whole 
without deforming it.) As an example, we assume that the droplet interface is pinned at ip\ (another example, force 
balance by suitable localized counterstresses, is studied in Ref. [36]): 



u(ip — ipi) = 0. 

The general solution of the inhomogeneous Legendre equation (B7) is 



u{iP) = AP{iP) + BQ{iP) + -nR 



R 



drop 



irop 



7 Jv 

where A and B are integration constants and 

P(ip) := cos^>, 

Q(ip) := 1 + cos ip In tan 



da sina [P(a)Q{ip) - P(iP)Q(a)]Il(a), 



(BIO) 



(Bll) 



ip 



(B12) 
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are solutions of the homogeneous Legendre equation. For notational simplicity, evaluation of a function at ipo(ipi) will 
be denoted by the subscript om- 
From Eq. (BIO) one obtains 



:M-Rdrop = -AP 1 -BQ 1 . 



From the other boundary condition (B8) it follows that 



£F-Rdrop — APq + BQ' Q + 



r>2 

^drop 

7 



da sina W (a)H(a) , 



where we have defined the auxiliary function 

W(a) := [P(a)Q' -P^Q(a)\. 
Finally in order to impose the integral constraint (B9) we employ the following identities: 

ci>i f>Pi 



I dipsinip / da sin a[P(a)Q(ip) - P(ip)Q(a)]U(a) = 

Jipo Jlp 

da sina U(a) / dtp smip[P(a)Q(i>) - P{ip)Q{a)] = 



(B13) 



(B14) 



(B15) 



(B16) 



/ 



da sin a II(cr) — 



W(a) 



Wo 



The first equality follows from interchanging the order of integration. The last equality follows most easily by noticing 
(i) that the ^-integral in the second line is, as a function of <r, a solution of the differential equation (B7) with u 
replaced by this integral, /i-Rdrop replaced by —1, II = 0, and the analogues of the boundary conditions uq = = u' , 
and (ii) that the function (W(a)/Wo — l)/2 is the solution to the same differential equation with the same boundary 
conditions. With these identities Eq. (B9) turns into 



= AV + BQ - - / da sina 11(a) 
where we have introduced the numerical constants 



W(a) 



W 



(B17) 



V :-- 
Q ■-- 



di) sinV> [P(ip) - Pi] 



/ di 

Jtflg 



dip simp [Q(ip) - Qi}. 



(B18) 



We note that the constants A and B are determined by the two linear Eqs. (B14) and (B17). 

These results permit us to study the intermediate asymptotic regime ro !rc f,?' -C -Rdrop- In the expressions derived 
above, we have to take the limit i?d 

rop * while keeping fixed ipi , v — Rdvop^ &nd ro ? ref = ■RdropV'O (so that 
ipo, ip -C 1). In this limit, the constants V and Q are finite, while P(ip) — 1 and Q(ip) — ln(r/i?dro P ) — In 2 + 1. The 
integral over Id in Eq. (B14) simplifies as follows (provided that the stress field H(ip) decays sufficiently fast): 



/ da sincr W(a)U(a) = ds sin 

Jibn Jr 



W 



-^drop 



R, 



drop 



^drop 

ds sll 



Rd 



rop 



n 



rop 



Rd 



rop 



(B19) 



where we have used IF(s/i?dro P ) — -Rdrop/^o- The integrals over Id in Eqs. (B17) and (Bll) can be simplified in the 
same manner. From Eq. (B14) we obtain 



B ~ r (e F - e n ), 



(B20) 
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and Eq. (B17) yields 



Q 

A ~ -— r (e F - en). 



(B21) 



One finds indeed that the Lagrange multiplier /x determined by Eq. (B13) is linear in ep — en- Finally, the general 
solution (Bll) simplifies to 



u~yl(l-Pi) + S In 



-Rdr 



- In 2 + 1 - 



1 

7 Jr 



ds s n 



R 



drop 



r 



(B22) 



which renders Eq. (24) with the constant 



C = 2 exp 



(1-Pi)| + Qi-1 



(B23) 



This constant is finite and non-vanishing provided ipi is not close to w. The case -01 = ir is pathological because 
u'(ir) = for a smooth profile, so that the boundary condition (BIO) would overdetermine the problem. 
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